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ABSTRACT 


GCP-MOD2  is  an  extremely  rapid  means  for  performing  frequency 
analyses  of  equivalent  circuits  of  cascaded  two-ports.  An  operator 
describes  his  circuit  to  GCP-MOD  2,  using  an  abbreviated  plain 
English  vocabulary,  and  specifies  measurements  to  be  made  at  an 
appropriate  terminal  plane  for  frequencies  of  interest.  GCP  com¬ 
municates  laconically  with  the  operator  as  it  accepts  the  circuit 
description  and  any  other  instructions.  It  then  tabulates  the  re¬ 
sults  of  the  measurements,  and  awaits  further  instructions,  such 
as  circuit  modifications  or  different  measurements. 

A  sa  typical  example,  a  microwave  filter  circuit,  made  up  of  lengths 
of  transmission  line  and  five  resonators,  of  both  lumped  and  dis¬ 
tributed  elements,  can  be  described  to  GCP  and  evaluated  at  a  dozen 
frequencies  for  such  measurements  as  insertion  loss,  transmission 
phase-shift,  and  input  impedance,  all  in  about  10  minutes. 
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Al'TOMATIt  MIC  ROWAVE  CIRCUIT  \NALYS1S  WITH  GCP- MOD  Z 


I  l)ES(  RIPTION  OF  GCP- MOD  2 

\.  Introduction 

(i(’P-\U)I)  2  performs  frequency  analysis  on  equivalent  circuits  of  cascaded  linear  two-ports. 
It  does  this  by  generating  the  transfer  matrices  of  two-ports  described  to  it,  and  multiplying 
them  to  arrive  at  an  overall  transfer  matrix  for  the  cascade.  The  elements  of  the  final  matrix 
are  manipulated  and  operated  on  to  determine  "measured"  quantities  at  each  frequency. 

To  work  with  GCP,  an  engineer  describes  his  circuit  using  abbreviated  plain- English  words 
given  in  the  SECTION  MENU.  He  directs  that  various  measurements  be  made  on  the  circuit, 
using  words  from  the  MEASUREMENT  MENU.  He  may  change  the  circuit,  tune  elements  in 
the  circ  uit,  set  up  additional  circuits,  call  for  measurements,  and  terminate  operation  hy  using 
appropriate  words  from  the  list  of  COMMANDS, 

The  philosophy  on  which  the  design  of  GCP  is  based  is  that  the  typical  user  is  a  research- 
oriented  engineer,  rather  than  a  full-time  microwave-circuit  designer,  and  that  he  may  use 
GCP  rather  infrequently,  but  intensely.  It  is  assumed  that  the  user  does  not  want  to  sacrifice 
much  time  in  learning  or  relearning  to  use  GCP.  On  the  other  hand,  he  does  not  want  to  sacri¬ 
fice  speed  or  flexibility  of  operation  for  a  large  amount  of  instructional  information  from  GCP 
itself. 

The  compromise  reached  is  a  vocabulary  of  words  of  up  to  four  letters,  which  are  either 
truncated  English-language  words  or  other  mnemonic  abbreviations.  The  four-letter  maximum 
was  chosen  on  the  basis  of  giving  the  most  information  for  the  least  redundancy,  without  using 
excessive  typewriter  printing  time. 

Another  consideration  is  the  relative  control  of  the  action;  is  GCP  controlling  the  course 
of  action,  or  is  the  user  in  command?  This  is  another  compromise  situation,  but  the  COMMAND 
LIST  gives  the  user  options  to  redirect  the  course  of  analysis  while  GCP  controls  the  details 
of  the  specific  action  chosen  by  the  user. 

13.  Describing  a  Circuit 

1.  Numbering  Branches  and  Sections 

The  circuit  to  he  described  must  he  capable  of  representation  as  a  cascade  of  two-ports. 
Allowable  two-ports,  termed  "sections"  are  given  in  the  SECTION  MENU. 

In  setting  up  a  cascade  of  sections,  it  is  necessary  to  observe  the  following  system  of 
numbering,  in  order  that  each  section  used  have  the  proper  location  with  respect  to  others. 
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Each  cascade  setup  is  termed  a  "branch,"  and  is  identified  by  a  two-digit  number:  01, 

02,  ....  09,  10.  A  maximum  of  10  separate  branches  is  permissible  at  present. 

Each  section  setup  is  identified  by  the  number  of  the  branch  of  which  it  is  a  part,  and 
another  two-digit  number  giving  its  position  in  the  cascade.  For  example,  01,  05  refers  to  the 
fifth  section  in  the  first  cascade.  A  maximum  of  50  sections  is  permissible  in  each  branch. 

In  consideration  of  the  convention  that  a 
generator  is  on  the  left,  and  a  load  is  on  the 
right,  the  numbering  proceeds  from  right  to 
left,  and  includes  the  load,  if  any.  Figure  1 
illustrates  the  numbering. 

The  terminals  on  the  left  of  each  section 
bear  the  same  designating  number  as  the 
section.  This  convention  is  used  in  specify¬ 
ing  the  port  at  which  a  measurement  is  to  be 
made. 

GCP  indicates  its  readiness  to  operate  by  printing  01,  01.  The  user  then  types  the  name 
of  the  first  section  of  the  first  branch,  preceding  the  name  by  SHUN  or  SERI  if  the  section  is 
made  of  a  two-terminal  element  connected  in  shunt  or  series  to  form  a  two-port.  The  names 
and  descriptions  of  sections  are  found  in  the  SECTION  MENU. 

The  example  below  shows  how  branch  01  of  the  preceding  illustration  is  described  to  GCP: 


0.565  INCHES  1 


1  0  9  nH 


Fig.  1.  Example  of  section  numbering. 


GCP 

EXECUTION  BEGINS 

01, 01 

RES!  (OHMS) 

01,  02 

INDU  (NH) 

01,  03 

ZO,  (LENG  (IN),  DIEL) 
01,  04 

RESI  (OHMS) 

01,  05 


USER 

shun,  resi 
100 

seri,  indu 
.9 
teml 

77.565,  2.00 
shun,  resi 
50 


No  voltage  or  current  sources  are  used  with  GCP-MOD  2,  because  it  is  limited  to  the 
analysis  of  linear  two-ports.  Driving  point  and  transfer  ratios  are  adequate  to  characterize  the 
performance  of  such  circuits.  For  this  reason,  it  is  equally  valid  to  number  sections  within  a 
branch  from  left  to  right  or  from  right  to  left  as  shown  in  the  illustration.  However,  care  must 
be  taken  in  evaluating  measurements  and  in  using  nonreciprocal  sections,  since  right  to  left 
numbering  has  been  assumed  in  developing  GCP. 


2.  Use  of  Operator  Sections 

At  the  end  of  the  SECTION  MENU  are  the  OPERATOR  SECTIONS.  These  are  special 
operations  occupying  one  section  of  a  branch.  First  given  are  SHUN  JUNC  and  SERI  JUNC 
which  connect  a  side  branch  (with  its  own  branch  number)  to  the  main  branch,  in  a  shunt  or 
series  connection,  as  desired.  It  is  to  be  noted  that  it  is  the  side  branch  end  having  the 


higher  section  number  that  is  connected  at  the  main  branch.  The  lower  section  number  end 
(p,  00)  must  always  remain  unconnected.  Thus,  side  branches  may  not  connect  to  other  side 
branches,  nor  return  to  the  main  branch.  That  is,  they  may  not  form  closed  loops  in  the  cir¬ 
cuit  topology.  However,  a  side  branch  may  have  one  or  more  subsidiary  branches,  but  these 
branches-off-branches  may  not  branch  further.  It  is  important  to  note  that  any  set  of  measure¬ 
ments  may  be  made  only  within  a  single  branch.  Thus,  branches  may  not  be  connected  end-for- 
end,  nor  may  transfer  measurements  be  made  between  a  generator  in  one  branch  and  a  load  in 
another. 

One  main  branch  may  have  many  side  branches.  If  identical  side  branches  are  needed,  a 
single  side  branch  may  be  made  up  and  specified  at  each  junction. 

\  special  pair  of  operator  sections  allows  a  side  branch  to  rejoin  the  main  branch,  so  that 
the  two  run  parallel.  This  is  illustrated  in  OPERATOR  SECTIONS.  A  parallel  branch  is  in¬ 
itiated  bv  the  command  PARE  p,  where  p  is  the  two-digit  branch  number  of  the  side  branch. 
Farther  along  the  main  branch,  the  side  branch  is  returned  by  using  a  section  called  COMB, 
for  combine.  Observe  that  the  lower-numbered  section  ends  of  the  two  branches  involved  are 
connected  together  in  parallel  at  PARB,  and  the  higher-numbered  ends  are  connected  in  parallel 
at  COMB.  Thus,  the  branches  not  only  run  parallel,  they  are  also  in  a  parallel  connection  at 
each  end. 

A  series  connection  is  possible,  however,  at  either  or  both  ends.  Section  Il-G-2  describes 
the  method  for  making  series  connections  of  parallel  branches. 

Another  operator  section  is  1MAG,  for  image,  which  allows  the  user  to  select  any  cascade 
of  sections  in  a  branch,  and  duplicate  it  in  reverse  farther  along  the  branch.  This  is  a  valuable 
capability  for  multi-element  filters,  which  are  often  symmetrical  about  the  center  of  the  struc¬ 
ture.  The  duplicated  part  of  the  cascade  occupies  only  one  section,  even  though  the  original 
occupied  a  number  of  sections. 


C\  Physical  Units 

As  set  up  at  present,  GCP- MOD  2  uses  the  following  set  of  units  because  they  are  convenient 
at  microwave  frequencies: 


Length 

Frequency 

Impedance 

Admittance 

Inductance 

Capacitance 


inches 

GHz  (109llz) 

ohms 

mhos 

nil  (10_9H) 
pK  (10‘12  F) 


I).  Making  Measurements 

When  the  operator  wishes  to  have  measurements  made  on  the  circuit  he  has  described  to 
GCP,  he  types  MEAS.  GCP  responds  by  printing  BR,  SE,  TYPE,  thus  asking  the  operator  for 
branch  and  section  numbers  of  the  terminal  plane  at  which  the  measurement  is  to  be  initiated, 
and  the  type  of  measurement,  given  by  the  name  in  the  MEASUREMENT  MENU.  The  operator 
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may  specify  up  to  three  types  of  measurements  to  be  made,  provided  all  are  initiated  at  the 
same  terminal  plane. 

There  are  three  basic  types  of  measurements.  First,  there  are  driving  point  measure¬ 
ments  of  impedance  or  admittance,  which  can  be  made  at  any  terminal  plane  looking  in  either 
direction.  Next,  there  are  driving  point  measurements  which  must  be  made  adjacent  to  a  purely 
resistive  termination  in  order  to  provide  a  reference  impedance  level.  This  includes  such 
measurements  as  VSWR  and  reflection  loss.  Finally,  there  are  transfer  measurements,  such 
as  transducer  loss,  or  transmission  phase  shift,  which  require  reference  impedance  levels  at 
both  ends  of  the  overall  cascade.  Thus,  transfer  measurements  are  made  only  on  circuits 
having  resistance  terminations  at  both  ends,  at  terminal  planes  adjacent  to  the  terminations.  It 
is  necessary  to  specify  only  the  terminal  plane  assumed  to  be  the  generator  port  for  that  meas¬ 
urement  GCP  assumes  the  other  termination  is  the  load  for  that  measurement. 

To  return  to  the  circuit  example,  assume  that  the  operator  wishes  to  measure  input  ad¬ 
mittance,  VSWR,  and  transducer  loss  for  branch  01  of  the  illustration,  assuming  the  generator 
port  to  be  01, 03.  The  GCP  user  interaction  is  as  follows: 

GCP  USER 


meas 


BR,  SE,  TYPE 


01,  03,  ylod,  vswr,  trio 


FREQ,  (MIN,  MAX,  1NCR)  4,  5,  0.2 


FREQ  (GHZ) 

ADMITTANCE 

VSWR 

TRAND  LOSS  (DB) 

4.0000 

0.0180 

0.0025 

1.1861 

0.0316 

4.2000 

0.0186 

0.0017 

1.1196 

0.0139 

4.4000 

0.0190 

0.0006 

1.0611 

0.0038 

4.6000 

0.0192 

-0.0005 

1.0501 

0.0026 

4.8000 

0.0191 

-0.0018 

1.1064 

0.0111 

5.0000 

0.0187 

-0.0030 

1.1808 

0.0300 

NEXT 


Measurement  of  individual  branches  may  be  made  at  any  time,  but  when  a  circuit  has  side 
branches  (using  JUNC)  or  has  parallel  branches  (using  PARB  and  COMB),  measurements  for 
the  entire  circuit  may  be  made  only  on  the  main  branch.  When  GCP  operates  on  the  main  branch 
it  "sees"  the  side  or  parallel  branches,  but  does  not  "see"  the  main  branch  when  operating  on 
a  side  or  parallel  branch. 

MEASUREMENT  DEFINITIONS  are  given  in  Sec.  II-C. 


E.  Use  of  Commands 

1.  Normal  Commands 

The  action  of  GCP-MOD  2  can  be  directed  by  the  operator  using  the  command  words  given 
in  the  COMMAND  LIST. 

El, MO  and  SEMO  allow  modification  of  an  element-value  and  a  section-type,  respectively. 
BRAN  n,  where  n  is  the  two-digit  cascade  number,  allows  the  operator  to  work  with  any 
particular  branch  called.  More  than  one  branch  might  be  needed  for  any  of  the  following  reasons 
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The  operator  may  wish  to  make  a  side  calculation,  or  to  try  out  a  subcircuit  before  incorporating 
it  in  the  main  branch,  or  the  main  branch  may  have  junctions  connecting  to  other  branches. 

The  command  MEAS  directs  GCP  to  accept  instructions  taken  from  the  MEASUREMENT 
MENU  described  previously. 

The  command  TUNE  is  used  when  the  operator  wishes  to  observe  the  changes  in  some 
measured  quantity  as  the  values  of  some  circuit  parameter  are  changed,  f  or  example,  the 
transducer  loss  of  a  filter  may  be  observed  at  a  critical  frequency  as  the  length  of  a  coupled 
stub  is  varied  over  some  range  of  values.  In  TUNE,  only  one  element  may  be  varied,  and  up 
to  three  types  of  measurements  may  be  made  at  only  one  frequency.  The  command  OFFL  is 
used  to  print  the  measured  data  on  an  offline  printer.  The  data  are  written  in  a  file  called 
GCP  DATA.  After  terminating  GCP  with  a  QUIT  command,  the  CMS  command  offline  print  gcp 
data  should  be  issued. 

OFFL  stays  in  effect  until  the  command  ONLI  is  issued.  ONL1  will  resume  the  typing  of 
measured  data  on  the  console. 

The  command  QUIT  is  used  to  terminate  operation  of  GCP.  The  information  or  data  of  the 
circuits  worked  with  are  not  erased  by  this  command.  Thus,  the  data  of  the  previous  operation 
are  available  when  GCP- MOD  2  is  turned  on  again.  To  reinstate  this  data,  it  is  necessary  to 
use  the  command  RESE,  for  reset,  when  GCP  prints  01,  01.  If  a  new  section  is  entered  when 
01,  01  is  printed,  the  file  of  previous  data  is  erased.  RESE  is  also  used  when  GCP  is  shut  down 
inadvertently  during  an  operating  session.  In  this  case,  it  should  be  noted  that  RESE  returns 
GCP  to  the  next  to  last  branch  previously  used,  rather  than  to  the  last  section.  This  is  done 
because  shutdown  may  have  occurred  after  the  last  section  had  been  entered,  but  before  it  had 
been  properly  filed.  Thus,  after  RESE,  it  is  necessary  to  re-enter  the  last  section  of  the 
branch  involved. 

When  GCP  has  completed  a  given  task,  it  prints  NEXT,  to  indicate  it  is  ready  to  accept  new 
instructions. 

2.  Diagnostic  Commands 

There  are  two  basic  diagnostic  commands  by  means  of  which  the  user  may  inspect  the  data 
on  which  GCP  is  operating.  The  first  of  these  is  LIST,  which  causes  GCP  to  print  out  the 
section  types  and  element  values  for  all  branches  and  sections,  or  only  specified  ones.  The 
numbers  printed  out  under  LIST  can  be  identified  by  referring  to  Sec.  II-D  under  "LIST”.  The 
other  diagnostic  is  MATR,  which  causes  GCP  to  print  out  the  values  of  the  transfer  matrix 
elements  (at  the  last  frequency  of  measurement)  of  all  or  some  of  the  sections  used  in  that 
sess  ion. 


K.  Options 

Different  areas  of  application  may  require  different  SECTION  MENUS  and  MEASUREMENT 
MENUS.  The  generation  of  special  menus  is  left  to  the  user. 

Anticipated  additions  to  GCP  will  provide  graphical  output,  and  an  optimization  capability. 

G.  Operating  Techniques 
1 .  Techniques 

A  typical  design  task  involves  evaluation  of  the  performance  of  a  proposed  circuit,  com¬ 
parison  of  the  performance  with  design  goals,  and  modification  of  the  circuit  in  hopes  of  improving 
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its  performance.  The  cycle  of  evaluation,  comparison,  modification,  and  re-evaluation  is  re¬ 
peated  until  a  realistic  circuit  having  satisfactory  performance  has  been  achieved,  or  found  to 
be  unobtainable. 

Thus,  in  the  interest  of  efficient  and  effective  use  of  GCP,  the  following  technique  is  offered 
(introduce  no  more  uncertainty  in  the  stated  problem  than  is  absolutely  necessary!): 

If  the  circuit  can  be  broken  down  into  subcircuits  which  can  be  optimized  separately,  then 
these  subcircuits  should  be  worked  on  separately. 

If  possible,  start  with  a  circuit  of  known  response  that  can  be  converted  into  the  actual 
circuit,  then  measure,  modify  and  optimize  one  element  at  a  time.  For  example,  consider  a 
filter,  bound  by  physical  constraints,  to  be  made  in  part  of  lumped  elements,  and  in  part  of 
distributed  elements.  Start  with  an  entirely  lumped  element  prototype  having  the  performance 
required,  since  lumped  element  filters  have  been  completely  determined  by  analytical  means, 
and  their  properties  published  widely.  Then,  investigate  the  change  in  performance  when  only 
one  lumped  element  section  is  replaced  by  the  required  distributed  section.  Optimize  perform¬ 
ance  for  this  single  change  before  going  on  to  the  next. 

With  some  problems,  it  is  possible  to  start  with  a  simple  but  meaningful  approximation  to 
the  actual  circuit,  using  EMPT  for  sections  which  will  later  be  required  to  incorporate  para- 
sitics.  Get  the  simple  approximation  to  work,  and  then  add  the  parasitics  one  at  a  time,  op¬ 
timizing  at  each  step. 

In  the  interest  of  simplicity  and  speed,  the  Operator  Command  IMAG  is  included  to  duplicate 
a  cascade  of  sections  for  a  symmetrical  network.  Observe  that  the  mirror  image  of  many 
sections  can  be  duplicated  using  only  one  section. 

Notice,  also,  that  if  a  number  of  identical  branches  is  needed  in  a  circuit,  it  is  only  nec¬ 
essary  to  define  that  branch  once,  and  then  specify  it  by  the  one  branch  number  at  each  junction 
(JUNC). 


2.  Series  Connection  for  a  Parallel  Branch 

The  pair  of  commands  PARB  and  COMB  permits  the  parallel  connection  of  a  branch  running 
parallel  to  the  main  branch.  The  connection  at  either  or  both  ends  can  be  made  a  series  con¬ 
nection  by  the  following  artifice.  Each  of  the  three  proper  circuit  sections,  which  would  normally 
connect  to  PARB  or  COMB  for  a  parallel  connection,  must  instead  be  required  to  be  in  cascade 
with  a  section-type  INTE.  This  makes  a  series  junction  because  INTE  interchanges  voltage  and 
current.  PARB  or  COMB,  which  in  effect  add  currents  to  make  a  parallel  connection,  will  then 
be  adding  what  are  called  currents,  but  which  are  actually  voltages,  thus  forming  a  series 
connection. 

3.  Transmission  Line  Loss 

GCP  does  not  have  a  section-type  for  transmission  lines  with  loss,  but  loss  in  transmission 
lines  or  waveguide  can  be  handled  accurately  by  the  following  method. 

If  we  take  the  transfer  matrix  for  lossy  transmission  line, 

coshyf  20  sinhyf 

-YOsinhyf  coshyf 
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where  y  <y  +  j ,  the  complex  propagation  factor,  and  4  is  the  length  of  line,  and  either  pre¬ 
multiply  or  post-multiply  by  the  inverse  matrix  of  a  lossless  transmission  line, 

cos/34  —  jZO  sin/34" 

.  -  j  YO  s  in  /34  c  os  £4 

which  means  that  we  are  removing  a  lossless  line  from  either  end,  we  are  left  with 
cosh  oi  It  ZOsinho4" 

.  YO  sinh  cv4  cosha4 


This  is  the  transfer  matrix  for  a  zero-length  attenuator  of  characteristic  impedance  ZO,  and 
loss  of  L  8,686  o?4  decibels.  For  a  typical  transmission  line  we  can  compute  or  measure  both 
ZO  and  o4,  or  L(dB).  Thus  a  lossy  line  can  be  treated  as  a  lossless  line  with  a  matched  atten¬ 
uator  at  one  end  or  the  other,  as  convenient. 

4.  Saving  a  GCP  Session  for  Future  Use 

If  GCP  is  inoperative  for  a  period  of  time,  overnight  for  example,  the  previous  data  are 
stored,  and  can  he  recalled  by  typing  RESE  in  response  to  GCP  typing  01,  01,  If,  however,  it 
is  desired  to  store  data  from  a  particular  session  without  disabling  GCP  for  other  use,  then  it 
is  necessary  to  ALTER  the  file  GCP  STORAGE  to  another  file  name  and  filetype,  or  have  it 
printed  out  on  cards  or  tape.  The  original  data  may  be  reinstated  with  GCP  by  subsequently 
renaming  the  second  file  back  to  GCP  STORAGE,  and  calling  RESE, 

5.  User-Generated  Sections 

Special  dummy  sections  (DUMA,  DUMB,  DUMC,  DUMD)  have  heen  incorporated  in  the 
SECTION  MENU  to  allow'  for  user-generated  special  sections.  These  are  difficult  to  implement 
w  ithout  a  thorough  knowledge  of  the  internal  details  of  GCP. 

II  MENUS 

A.  Section  MENU 

1  Circuit  Sections  (n  is  the  two-digit  branch  number,  and  k  is  the  two-digit 
number  of  the  section  to  be  introduced) 


ELEMENTARY  SECTIONS 


GCP  NAME 
AND  SYMBOL 


GCP 

WILL  TYPE 

REQUEST  OF  SECTION  TRANSFER  MATRIX 


EMPT 

o - o 


o - o 

n,K  n,k-l 


Straight  through 
connection  oc¬ 
cupying  one 
section 


1  0 
0  1 
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ELEMENTARY  SECTIONS  (Continued) 


CCP  NAME 
AND  SYMBOL 


GCP 

WILL  TYPE 

REQUEST  OK  SECTION  TRANSFER  MATRIX 


shor 


o 

o 


Short  circuits 
load  end  of 
cascade 


SIIOR  -  Short  circuits  0  1 

generator  end 

o  of  cascade  (  Ll  0. 

- o 

n,k-l 


INTE 


n,k  n,k-l 


Interchanges 
voltage  and 
current 


Inverts  polarity 
of  terminals 


-  1 

0 


0 

-1 


n,k  n,k-l 


XFMR 


TURNS  RATIO  Ideal 

transformer 


N  0 
0  1/N. 


INVE 


ZO(OHMS)  Ideal  imit- 

tance  inverter 
with  90°  phase 
shift 


0  jZO 
jYO  0  . 


ATTN 


ZO,  Ideal  attenuator 

LOSS(I)B) 


U'-’f1  77  > 

vn  (L  j)  j  u.  ‘  |  ) 

J  _  10LOSS(IJm/20 
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ELEMENTARY  SECTIONS  (Continued) 


GCP  NAME 
ANI)  SYMBOL 


GCP 

WILL 

TYPE 

REQUEST 

OF  SECTION 

TRANSFER  MATRIX 

ISOL 


ZO,  LOSS, 
REVL(DB) 


Ideal  isolator 


ill 


¥  H. I- 


i  .IT » 


1. 1? 


1  )  y  (U 


I.Ii 


,  tr  ,  AI  ()SS(l)Hl/20  .  „  .  rtU.E  VI .  (Dili  20 

LF  10  I.R  10 


SECTIONS  MADE  OF  TWO-TERMINAL  CIRCUITS 


SHUN  NAME  where  NAME  Info  needed  Admittance  in 

stands  for  GCP  name  of  for  circuit  SHUNT 

two-terminal  circuit  con-  named 
nected  in  shunt  to  forma¬ 
tion  n,  k 

o - f - o 


2-TERMINAL 

CIRCUIT 


O - * - O 

n,  k  n,k-1 


1  0 
Y  1 


SERI  NAME  where  NAME  Info  needed  Impedance  in 

stands  for  GCP  name  of  for  circuit  SERIES 

two-terminal  circuit  con-  named 
nected  in  series  to  form 
section  n,  k 


2-TERMINAL 

CIRCUIT 

1  Z 
0  1 


o - o 

n,k  n,k-1 


LUMP  CIRCUIT  ELEMENT  SEC  TIONS 
(SHUN  or  SERI  must  precede  NAME  of  two-terminal  circuits) 


GCP 

GCP  NAME 

WILL 

TYPE 

AND  SYMBOL 

REQUEST 

OF  SECTION 

IMMITTANCE 

RESI 

RES  (OHMS) 

Resistance 

Z  =  R 
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LUMP  CIRCUIT  ELEMENT  SECTIONS  (Continued) 

(SHI  N  or  SERI  must  precede  NAME  of  two-terminal  circuits) 


GCP  NAME 
ANL)  SYMBOL 


GCP 

WILL  TYPE 

REQUEST  OF  SECTION  IMMITTANCE 


INDU 


IND  (NH) 


Inductance  Z  -  jwL 


CAPA  CAP  (PF)  Capacitance  Y  =  jcoC 

i 

T 


SRLC 


SRES 


RES  (OHMS),  Series  resonant  Z  =  R  +  j(wL  —  — tt  ) 
IND  (NH), 

CAP  (PF) 


FO(GHZ),  X,  Series  resonant  „  X  .  F  FO, 
1/Q  Z  Q  +  jX(FO  “ 


where 


FO  =  Resonant  freq. 
x  27:  FOL,  slope  parameter 
Q  =  2tt  FOL/R 


o 

PR  LC 


PRES 

6 


RES  (OHMS),  Parallel 

IND  (NH),  resonant 

CAP  (PF) 

FO(GHZ),  BP,  Parallel 

1/Q  resonant 


I  ^  Is 

R  +J<wC"^L) 


Y  =  +  jBP(—  —  -EQ\ 

Y  Q  +  FQ  F 


where 


FO  =  Resonant  freq. 

BP  27r  FOC,  slope  parameter 
Q  -  2tt  FOCR 
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TRANSMISSION  LINE  TWO- PORT  SECTIONS 
(SHUN  or  SERI  must  precede  NAME  of  two-terminal  circuits) 


GCP  NAME 
AND  SYMBOL 


GCP 

WILL 

TYPE 

REQUEST 

OF  SECTION 

IMMITTANCE 

TEMS 

u 


Z O,  LENG 
(IN.),  DIEL 


Short-circuited 
TEM  transmis¬ 
sion  line 


Z  *  / < ) tan 


2jrF<I.KN(!)  (1)IK1.»*  2 
c 


DIKI.  =  Itclativo  dielectric  constant 


TEML 


ZO,  LENG  Open  circuited  _  t (Dim.)1  2 

(IN.),  DIEL  TEM  transmis-  * J 

sion  line 


GCP  NAME 
AND  SYMBOL 


TEML 


n,k  n,k-1 


GCP 

WILL  TYPE 

REQUEST  OF  SECTION  TRANSFER  MATRIX 


ZO,  LENG 
(IN.),  DIEL 


TEM  transmis-  cos  0  jZOsinO 

sion  line 
(lossless ) 

.  j  Y o  .sinO  cos  6 


where 


G 


2jtF(LENG)  (DIED1  2 

c 


and 


c  11.8  inches  nsec 


H  \DL 


-H" 

n,  k 

n, k  - 

J)LOD,  DGEN, 
B,  DIEL 

DLOD  DIA, 
Nearest  Load, 

DGEN  DIA, 
Nearest  GEN, 

B  Spacing  of 
Plates 


See  Appendix. 
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WAVEGUIDE  ELEMENT  SECTIONS* 


GCP  NAME 
AND  SYMBOL 


GCP 

WILL  TYPE 

REQUEST  OF  SECTION  TRANSFER  MATRIX 


TEWG 


WIDE,  NARO, 
LENG,  DIEL 


TE-mode  rec-  p  cos  0  jZO  sinG 

tangular 

waveguide 

_  jYO  sin 0  cos  0 


o 


i  ( 


1 1  ,K\( ;  1 1 1  >1  r:i , i 


2  (WIDE)  ( 1 ) 1 1; i . » 1  2 

....  24  OrlNAKOt  I-  _ _ 

/u  fY*  2  *  2  1  '2 

(WIDE)  (l)IKI  i  (I  1(1 


WIDE,  NARO,  Waveguide 
DIEL  termination 


'  1  0 

YO  1 


YO  l/ZO  of  TEWG. 


Urn  |  TEWS 


WIDE,  NARO,  Short-  Z  =  jZO  tan  0 

LENG,  DIEL  circuited 

TE1 0 
waveguide 


*  The  definition  of  characteristic  impedance  in  column  4  is  correct  for  coupling  to  a  probet 
or  packaged  diode t  mounted  across  the  narrow  dimension.  It  is  also  suitable  for  cascaded 
waveguides  having  different  dielectrics  or  different  narrow  dimensions  (except  for  fringing 
capacitance ).  § 

For  cascaded  waveguides  of  differing  wide  dimensions,  both  an  impedance  transformation  and 
fringing  susceptance  must  be  accounted  for.JT 

tS.  A.  Schelkunoff,  Electromagnetic  Waves  (D.  Van  Nostrand,  Inc.,  Princeton,  New  .Jersey, 

1943),  Sec.  12.4,  pp.  494-496. 

JW.  J.  Getsinger,  "The  Packaged  and  Mounted  Diode  as  a  Microwave  Circuit,"  IEEE  Trans,  on 
Microwave  Theory  and  Techniques,  MTT-14,  58-69  (February  1966). 

§N.  Marcuvitz,  Waveguide  Handbook,  Rad.  Lab.  Series  (McGraw-Hill,  New  York,  1951),  Vol.10, 
Sec.  5.26,  pp.  307-310. 

5 Ibid.,  Sec.  5.24,  pp.  296-302. 
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WAVEGUIDE  ELEMENT  SECTIONS  (Continued) 


GCP  NAME 
AND  SYMBOL 


GCP 

WILL 

TYPE 

REQUEST 

OF  SECTION 

TRANSFER  MATRIX 

WISO 


n,k  -  n,k-1 


WIDE,  NAHO,  Waveguide 

DIEL,  LOSS,  isolator 

HEVL(DB) 


i  ....  1 


1 M  1 


VO  (M _ 

2  11,1  IK 


1  1  . 
i  -  d.i* 


U  »  10'  OSSU.HI/20  ,  1{  t0 


n< 1 

HK\  1 .  1 1  )H>  20 


The  GCP  names  for  user-initiated  sections  are: 

DUMA  DUMC 

DUMB  DUMD 


2.  Operator  Sections 


SYMBOL 


o 


o- 

n,  k 


■O 


p,  00 


- O 

n,k-1 


GCP  WILL 
GCP  NAME  REQUEST 

SHUN  JUNC  BR,  SE 


TYPE  OF  SECTION 

SHUNT  connected  JUNCTION 
to  another  branch 


TL 

BRANCH 

P 


O  O 

n,k  n,k-1 


SERI  JUNC  BR,  SE 


SERIES  connected  JUNCTION 
to  another  branch 


O  On,k-l 


r~n 

I 


L T-pJ 

» n,  m  -t 


6  6  m 


PARB  p 


COMB 


Junction  nearer  the  load  for 
initiating  a  PARALLEL 
BRANCH  (numbered  p)  in 
PARALLEL  connection  with 
main  branch  (numbered  n) 

JUNCTION  nearer  the  gen¬ 
erator  for  ReCOMBINING 
parallel  branch  (numbered  p) 
in  parallel  connection  with  the 
main  branch  (numbered  n). 

Note*  m  —  1  >  k 

Series  connected  junctions  are 
handled  as  described  in  Sec.  I-G-2. 
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SYMBOL 

GCP  NAME 

GCP  WILL 
REQUEST 

TYPE  OF  SECTION 

-cnF*oFia- 

IMAG 
n,  d,  n,  i 

START,  END 
(BR,  SE) 

The  section  placed  at  n,  k  will 
be  the  mirror  image  of  the 
cascade  of  sections  between 
n,  d  and  n,  i. 

k,  k  +  1  n,k  n,  i,  n,  d 

k  >  i  >  d 

B.  Measurement  MENU 

These  measurements  may  be  called  only  after  giving  the  command  MEAS  (n  and  k  are 
the  two-digit  branch  and  section  numbers,  respectively,  of  the  terminal  plane  at  which  the 
measurement  is  to  be  initiated). 


1  Driving  Point  Measurements  —  Absolute 
These  measurements  may  be  made  at  any  terminal  plane. 


n,  k  ZLOD 


n,  k  ZGEN 


n,  k  YLOD 


n,  k  YGEN 


Impedance  at  terminal  plane  n,  k  looking 
toward  decreasing  section  numbers'. 

Impedance  at  terminal  plane  n,  k  looking 
toward  increasing  section  numbers 

Admittance  at  terminal  plane  n,  k  looking 
toward  decreasing  section  numbers. 

Admittance  at  terminal  plane  n,  k  looking 
toward  increasing  section  numbers. 


2.  Driving  Point  Measurements  —  Ratios 

These  measurements  must  be  made  at  a  terminal  plane  between  a  pure  resistance  termi¬ 
nation  and  the  remainder  of  the  circuit. 


n,  k  REFL 

n,  k  REFM 
n,  k  RFPH 
n,  k  VSWR 
n,  k  RFLO 
n,  k  RTLO 


Complex  reflection  coefficient  (Sll  or 
S22,  if  properly  terminated) 

Magnitude  of  reflection  coefficient 

Phase  of  reflection  coefficient  in  degrees 

Voltage  standing  wave  ratio 

Reflection  loss  in  decibels 

Return  loss  in  decibels 


3.  Transfer  Measurements 

These  measurements  require  resistance  terminations  at  both  ends  of  the  branch  being 
analyzed.  Measurements  must  be  made  at  a  terminal  plane  between  a  pure  resistance  termi¬ 
nation  and  the  remainder  of  the  circuit.  The  terminal  plane  at  which  the  measurement  is  called 
represents  the  generator  end  for  that  measurement. 

n,  k  TRAN  Complex  voltage  transmission  coefficient 

[SI  2  if  measured  at  load  end  (n,  01),  or 
S21  if  measured  at  generator  end  (n,  k); 
k  >  01,  if  properly  terminated  at  both  ends]. 
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n,  k  TRMG 
n,  k  TRPH 
n,  k  TRLO 
n,  k  SCAT 


Voltage  transmission  coefficient  magnitude 

Voltage  transmission  coefficient  phase 

Transducer  loss  in  decibels 

Scattering  matrix  elements  (no  other  meas¬ 
urement  may  be  made  at  the  same  time) 


C.  Measurement  Definitions 


n,  k  ZLOD 


AL 

CL 


n,  k  YLOD 


CL 

AL 


n,  k  ZGEN  = 


DG 

CG 


n,  k  YGEN  = 


CG 

DG 


n,  k  REEL  = 


AL  _  DG 
CL  CG 

A I  DG 
CL  CG 


AL  DG 

n-01REFL  -gr  I  jg 

CL  CG 


n,  k  REFM  -  REEL! 
n,  k  RFPH  2  arctan 


Im  REEL 


Re  REE 


?L  +  J 


(Re  REE)2  +  (Im  REEL)2 


n,  k  VSWR 


1  ±  REEM 

1  -  REFM 


n.kRELO  10  1og10 


1  -  REEM 


n,  k  RTLO  =  20  log4  Q 


n,  k  TRAN  = 


Si  2 
S21 


if  k  04 
if  k  >  01 


n,  k  TRAN 


2[§i  ‘  YLOD(n'01^  2 


AL  + 


DG 

CG 


CL 


2  [RG(n,  k  +  l)/RL(nt  Ol)]1/2 

A(n,  k)  +  RG(n,  k  +  4)  C(n,  k) 


n,  01  TRAN  - 


2  [RL(n,  01  )/RG(n,  k  +  l)]1  2  [A(n,  k)  D(n,  k)  -  B(n,  k)  C(n,  k)] 

D(n,  k)  4-  RL(n,  01)  C(n,  k) 


n,  k  TRMG  =  n,  k  |TRAN| 
n,  k  TRMG  -  n,  01  |TRAN| 
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n,  k  TRPH  =  2  arctan 


Im  TRAN 


Re  TRAN  + 


J{ Re  TRAN)2  +  (Im  TRAN)2 


n,  k  TRLO  =  20  log1Q 


1 

TRMG(n,  k) 


n,  01  TRLO  =  20  l°g1  0  TRMG(n,  01) 

n,  k  SCAT,  Sll  -  n,  k  REFL 
S21  -  n,  k  TRAN 

D.  Commands 

NORMAL 

COMMAND  GCP  WILL  REPLY 

GCP-MOD  2  Execution  begins 

01,  01 

SEMO  BR,  SE,  NEW  SECT. 

ELMO  BR,  SE,  ELEM 

NEW  VALUE 
BRAN  n  n,  k 

MEAS  BR,  SE,  TYPE 

FREQ(GHZ),  MIN,  MAX, 
INCR 

TUNE  BR,  SE,  ELEM 


SI  2  =  n,  01  TRAN 
S22  =  n,  01  REFL 


COMMANDS 

FUNCTION  -  OPERATOR 
RESPONSE 


Initiates  operation  of  GCP-MOD  2.  Type 
in  name  of  first  section  in  first  branch 
(e.g.,  SHUN  RESI)  or  RESE  (see  below) 
if  previous  work  is  to  be  recalled. 

SECTION  MODIFICATION.  Type  in 
branch  and  section  numbers,  and  name 
of  new  section  wanted. 

ELEMENT  MODIFICATION.  Type  in 
branch  and  section  numbers,  and  name 
of  parameter  to  be  changed,  e.g.,  01, 

05,  ZO. 

Type  in  new  value  wanted. 

Sends  GCP  to  Branch  indicated  by  n, 
the  two-digit  number  of  the  branch, 
k  is  the  two-digit  number  of  the  first 
section  not  previously  used  in  that 
branch.  Type  in  name  of  section  wanted 
for  n,  k,  or  command  for  action  on  that 
branch. 

Initiates  MEASUREMENT.  Type  in 
branch  and  section  numbers  of  terminal 
plane  where  measurement  is  to  be  made 
(terminal  plane  adjacent  to  assumed 
generator  for  transfer  measurements) 
and  names  (up  to  three)  of  measurements 
to  be  made. 

Type  in  frequency  range  and  increment. 


TUNES  (increments)  value  of  one  sec¬ 
tion  parameter  and  makes  single  fre¬ 
quency  performance  measurements  at 
each  increment.  Type  in  branch  and 
section  numbers  and  name  of  element 
to  be  tuned. 
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NORMAL  COMMANDS  (Continued) 


COMMAND  GCP  WILL  REPLY 

TUNE  BR,  SE,  MEAS 

ELEM(MIN,  MAX,  INCR.), 
FREQ(GHZ) 

OF  PL  NEXT 

ONLI  NEXT 

RESE  n,  k 

QUIT 


FUNCTION  -  OPERATOR 
RESPONSE 


Type  in  branch  and  section  numbers  of 
terminal  plane  and  names  (up  to  three) 
of  measurements  to  be  made. 

Type  in  range  and  increment  for  varia¬ 
tion  of  element  value,  and  frequency  at 
which  measurement  is  wanted. 

Prints  measured  data  in  a  file  named 
GCP  DATA.  The  CMS  command  offline 
print  gcp  data  must  be  issued  after 
terminating  GCP-MOD  2. 

This  command  remains  in  effect  until 
ONLI  is  issued. 

Resume  printing  of  measured  data  on 
the  console. 

RESETS  GCP  at  next  to  last  developed 
section.  Usually  used  to  reinstate  data 
file  after  intended  or  unintended  inter¬ 
ruption  in  operation  of  GCP.  Retype 
name  of  last  section,  or  type  command 
for  other  action  wanted.  Data  file  will 
be  erased  if  RESE  is  not  used  in  re¬ 
sponse  to  GCP  printing  01,  01  at 
initiation  of  operation  of  GCP. 

Terminates  operation  of  GCP,  but  does 
not  erase  data  file. 


COMMAND 


DIAGNOSTIC  COMMANDS" 
GCP  WILL  REQUEST 


FUNCTION 


LIST  BR,  or  BR,  SE 

LIST  ALL 


MATR  BR,  or  BR,  SE 

MATR  ALL 


Prints  out  section  name  and  element 
values  for  BR,  SE  specified,  or  if  only 
BR  is  specified,  prints  out  section 
names  and  element  values  for  all  sec¬ 
tions  in  that  branch,  or  if  ALL  is  used 
in  the  initial  command,  prints  out  all 
section  names  and  element  values  of  all 
branches  used  in  that  session. 

Element  values  can  be  identified  as  to 
parameter  involved  by  comparing  order 
on  print-out  to  order  given  below  under 
LIST. 

Prints  out  elements  of  transfer  matrix, 
for  the  last  frequency  measured,  for 
BR,  SE  specified;  or  if  only  BR  is 
specified,  prints  out  all  elements  of  the 
section  transfer  matrices  for  that  branch, 
or,  if  ALL  is  used  in  the  initial  com¬ 
mand,  prints  out  elements  of  all  trans¬ 
fer  matrices  of  all  branches  used  in  that 
session. 


*  These  commands  may  be  used  at  any  point  in  a  session  with  GCP  to  exhibit  the  data  GCP  has 
been  working  with. 
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LIST 


SERI,  REST 
OHMS 

SERI,  CAPA 
FARADS 


SHUN,  RESI 
MHOS 

SHUN,  CAPA 
FARADS 


TEML 

OHMS 

LENG(IN.) 

1/2 

(DIED  7 


SERI,  INDU 
HENRIES 

SERI,  TEMS 
OHMS 
LENG(IN-) 
(DIEL)1/2 

SERI,  PR LC 
MHOS 
FARADS 
HENRIES 

SERI,  SRLC 
OHMS 
FARADS 
HENRIES 

SERI,  PRES 
2ttFO 
B(MHOS) 
1/Q 

SERI,  SRES 
2ttFO 
X(OHMS) 
1/Q 

SERI,  TEML 
OHMS 
LENG(IN.) 
(DIED1  /  2 

SERI,  JUNG 
BR 
SE 
2 

SERI,  TEWS 
WIDE 
LENG(IN.) 
(DIEL)1/2 
NARO 


SHUN,  INDU 
HENRIES 

SHUN,  TEMS 
OHMS 
LENG(IN.) 
(DIEL)1/2 

SHUN,  PRLC 
MHOS 
FARADS 
HENRIES 

SHUN,  SRLC 
OHMS 
FARAD 
HENRIES 

SHUN,  PRES 
2ttFO 
B(MHOS) 
1/Q 

SHUN,  SRES 
2ttFO 
X(OHMS) 
i/Q 

SHUN,  TEML 
OHMS 
LENG(IN.) 
(DIEL)1/2 

SHUN,  JUNC 
BR 
SE 
1 

SHUN  TEWS 
WIDE 
LENG 
(DIEL)1/2 
NARO 


I  SOL 
OHMS 
FORL  -  DB 
RE  VI  DB 

ATTN 

OHMS 

DB 

XFMR 

TURNS  RATIO 

INVE 

OHMS 

RADL 
DGEN/2 
(DIEL)1/2 
B  -  B 
DLOD/2 

IMAG 

} START 

SE  ^  END 

PARB 

BR 

WISO 

WIDE 

FORL  DB 
REVL  DB 
(DIEL)1/2 
NARO 


WLOD 

WIDE 

0 

(DIEL)1/2 

NARO 
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III.  PH OG HAM  DESCRIPTION 


A  Introduction 

GCP-MOD2  is  a  CP/CMS  Fortran  program  designed  to  allow  users  without  programming 
knowledge  to  perform  circuit  analysis  on  a  wide  range  of  circuits.  The  following  description 
is  somewhat  of  a  verbal  flow  chart  in  which  the  flow  of  information  supplied  by  the  user  is  proc¬ 
essed  by  GCP. 

GCP  has  one  general  control  program  to  which  different  circuit  libraries  can  be  attached. 
The  circuit  libraries  are  text  libraries,  containing  subroutines  written  in  a  specified  form, 
which  convert  information  about  the  circuit  supplied  by  the  user  into  information  processed  by 
GCP. 


~7 16- 4-12613] 


04,01  — JQ5’03  05,02  05,01 


SERI 

JUNC 


04,02  “3HH 


-  h  -  h  -  -  02,01 


04,03  -  ~  -f  - 


04,04 


3 eh; 


SHUN 

JUNC 


-  | - ' - 02,02 
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Fig.  2.  Typical  circuit  layout. 

Figure  2  shows  schematically  the  way  the  circuits  arc  laid  out.  The  layout  consists  of  up 
to  ten  branches  with  each  branch  capable  of  handling  fifty  sections.  Each  section  contains  a 
circuit  element  or  elements,  when  more  than  one  element  can  be  conveniently  represented  in 
one  section.  Each  section  is  converted  into  a  transfer  matrix.  The  sections  are  identified  by 
branch  and  section  number  (for  example,  03,  05  refers  to  the  third  branch,  and  the  fifth  section). 

The  control  program.  Fig.  3,  consists  of  several  executive  modules  which  control  the  flow 
of  information.  These  executive  modules  function  automatically  depending  on  the  information 
given  by  the  user.  Information  given  by  the  user  is  stored  in  two  arrays.  One  array,  which 
contains  the  transfer  matrices,  is  ABCD  (10,  50,  4)  and  is  dimensioned  as  shown,  representing 
10  branches,  50  sections,  and  4  matrix  elements.  The  array  is  complex.  The  second  array, 
STOR  (10,  50,  4,  2,  2),  stores  information  pertaining  to  branch  number,  section  number,  element 
values,  section  type,  and  element  type,  respectively. 

The  STOR  array  is  also  written  into  the  data  file  with  data  set  reference  number  4  for  a 
permanent  record  of  the  information  in  case  of  system  failure  or  fatal  error  on  the  user’s  part. 
The  file  name  is  GCP  STORAGE. 
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Fig.  3.  Simplified  GCP  flow  diagram. 

These  two  arrays  are  in  a  labeled  common  named  INDEX.  There  are  several  labeled  com¬ 
mons  which  are  needed  to  transfer  information  between  subroutines. 

An  exec  program  is  used  to  start  GCP.  It  does  a  GLOBAL  LOADER  and  a  LOAD  GCP 
(NOMAP  XEQ),  loading  the  appropriate  text  libraries  and  GCP.  A  LOAD  MODULE  is  also 
available,  and  the  program  may  be  started  by  issuing  GCP-MOD2.  Once  the  program  is  started, 
all  arrays  are  initialized  to  zero.  This  feature  is  used  to  test  the  arrays  to  find  out  how  much 
of  each  is  used.  After  initializing,  the  first  executive  module  (section  exec,  Fig.  3)  is  entered. 
This  module  sets  the  branch  and  section  number  and  checks  to  make  sure  the  branch  and  section 
limits  are  not  exceeded,  then  writes  the  data  from  the  last  section  formed  into  GCP  storage  on 
the  disk  file,  using  a  direct  access  write.  The  section  exec  then  prints  the  next  available  branch 
and  section  number  on  the  typewritten  console.  Present  branch  and  section  numbers  are  stored 
in  the  integer  variables  N  and  M,  respectively. 
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The  user  types  in  a  section  type  or  a  command  and  the  information  is  read  into  two  integer 
variables,  TYPE  and  ELEM1  under  an  A4  format.  The  information  supplied  by  the  user  is  first 
compared  to  a  list  of  commands  using  logical  if  statements,  if  a  command  was  given,  the  ap¬ 
propriate  logical  if  executes  the  command.  If  a  section  request  is  given,  a  subroutine  SUB  is 
called.  This  subroutine  passes  the  arguments  TYPE  and  ELEM4.  The  subroutine  contains  a 
list  of  section  types,  and  the  request  is  compared  to  the  list  of  section  types  using  logical  lfs. 

If  no  match  is  found,  an  error  is  assumed  and  an  error  message  is  typed  on  the  typewritten 
console.  The  branch  and  section  numbers  are  reset  and  control  is  returned  to  the  section  excc 
where  the  branch  and  section  number  is  typed  out.  The  section  exec  is  then  ready  for  more 
input  information.  If  a  match  is  found,  the  appropriate  section  subroutine  is  called. 

B.  Section  Subroutines 

Section  subroutines  are  generally  set  up  in  the  following  manner 

The  subroutine  name  is  the  name  of  the  section  with  a  prefixed  Q  (i.e.,  TEML  is  QTEMI  ). 
Generally  no  arguments  arc  passed  except  one  or  more  return  statement  numbers  using  *  as 
defined  in  Fortran  IV. 

Complex  statements  for  the  various  complex  variables  are  required,  along  with  the  appro¬ 
priate  common  statements.  Labeled  common  INDEX  is  always  required  in  all  subroutines  since 
it  contains  N,  M  and  the  STOR  array.  These  variables  and  array  are  needed  to  identify  the 
section.  The  labeled  common  ZMEAS  is  generally  required  because  it  contains  section  and 
branch  information  needed  in  the  entry  section. 

Generally  a  NAMELIST  is  required  since  each  section  requires  specific  information  to  be 
supplied  by  the  user.  This  is  accomplished  by  a  write  statement  which  asks  for  the  information 
and  a  NAMELIST  READ  which  reads  it. 

Section  information  (length  of  a  section  of  transmission  line  for  example)  is  placed  in  the 
STOR  array.  Eight  storage  locations  are  available  for  this  data.  They  are  the  third  and  fourth 
dimensions  of  the  array. 

Finally  each  section  contains  an  identification  code  stored  in  STOR  (N,  M,  4,  2.  1).  This 
code  is  used  to  recall  the  subroutine. 

To  this  point,  the  subroutine  has  stored  the  information  pertaining  to  a  specific  section, 
and  control  passes  back  to  the  section  exec. 

Each  section  subroutine  contains  an  entry  point  The  entry  name  is  the  same  as  the  sub¬ 
routine  name,  but  with  a  Z  prefixed  instead  of  Q.  (i.e.,  TEML  would  be  ZTEML).  The  entry 
is  called  for  a  measurement  and  generates  the  transfer  matrix  using  information  stored  in  the 
STOR  array. 

C.  Commands 

All  of  the  remaining  subroutines  to  be  descrihed  pertain  to  commands.  The  commands  are 
few  in  number  and  are:  BRAN  XX,  RESE,  SEMO,  ELMO,  MEAS,  TUNE,  LIST,  MATR,  and 
QUIT. 

BRAN  XX,  where  XX  is  a  branch  numher,  either  starts  a  new  branch  if  the  branch  is  empty 
or  enters  the  branch  at  the  next  availahlc  section.  BRAN  and  the  branch  number  (XX)  are  read 
into  integer  variables  TYPE  (or  ACTION)  and  ELEM4  under  an  A4  format. 
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The  list  of  commands  is  scanned  using  logical  ifs  and  the  appropriate  if  statement  calls  tin* 
subroutine  QJLNC.  QJt  NC  tests  ELEM1  to  be  sure  \  branch  numl)(*r  was  given  (be.,  tests  to 
see  if  it  is  not  blank).  If  it  is  blank  an  error  message  is  typed  out  on  the  console  and  control 
is  passed  back  to  the  user  for  appropriate  action. 

If  no  error  is  detected,  ELEM1  must  be  changed  from  an  A  formatted  variable  to  an  1 
formatted  variable.  This  is  accompl ished  by  writing  it  into  a  file  with  data  reference  set  number 
1  and  reading  it  back  as  an  I  format. 

Esing  a  do  loop,  STOR  (El. EMI,  1,  1,  2,  1)  is  first  checked  to  see  if  it  is  zero  (this  location 
contains  the  program  identification  code).  If  it  is  zero,  this  indicates  that  it  is  an  empty  branch 
and  therefore  a  new  branch  is  being  called.  The  section  variable  M  is  set  to  zero  and  control 
retui  ned  to  the  section  exec  and  the  user. 

If  STOR  (ELEM1,  1,  1,  2,  1)  is  not  zero,  the  do  loop  continues  and  cycles  the  second  dimen¬ 
sion  until  a  zero  is  found  indicating  that  this  section  is  empty  (the  second  dimension  is  the  storage 
location  containing  section  numbers).  The  section  variable  is  set  to  this  number  and  control  is 
returned  to  the  section  exec,  whereupon  the  section  exec  types  the  branch  and  section  numbers 
at  the  console,  and  control  passes  to  the  user. 

The  command  RESE  allows  the  user  to  reset  GCP.  This  command  is  usually  given  at  the 
start  of  the  program  to  recall  all  information  into  active  storage.  The  command  was  imple¬ 
mented  to  allow  the  user  to  recall  previously  stored  information  in  case  of  a  system  crash  or 
fatal  user  error. 

As  previously  stated,  all  information  is  written  into  file  GCP  STORAGE.  The  RESE  com¬ 
mand  sets  the  file  pointci  to  the  top  of  the  file  and  sequentially  reads  this  direct  access  file  into 
the  STOR  array.  Since  the  information  is  read  originally  into  GCP  STORAGE,  one  section  behind 
the  present  section,  RESE  resets  to  that  section.  This  command  must  be  given  as  the  first 
command  if  a  reset  is  desired.  If  it  is  not  given,  the  GCP  storage  file  is  erased  using  CP/CMS 
CALL  ERASE  statement. 

RESE  may  be  given  anytime  during  execution  and  the  active  storage  will  be  reset  to  the 
previous  section 

The  SEMO  command  calls  a  subroutine  which  allows  the  user  to  modify  or  change  a  section. 
The  subroutine  asks  for  branch,  section,  and  section  type,  then  does  a  read  similar  to  the 
section  exec  read  for  a  section  request. 

STOR  (N,  M,  1.  2,  1)  is  checked  for  zero  to  see  if  a  section  exists.  If  STOR  (N,  M,  1,  2,  1) 
is  zero,  an  error  message  is  printed  and  control  passed  back  to  the  command  exec  and  the  user. 

If  STOR  (N,  M,  1,  2,  1)  is  not  zero,  a  list  of  section  types  is  scanned  using  logical  ifs.  If  a 
match  is  found,  the  appropriate  subroutine  is  called  and  the  STOR  array  updated  for  that  section. 

In  the  event  no  match  is  found,  an  error  message  is  printed  and  control  is  returned  to  the* 
command  exec  and  the  user  for  corrective  action 

If  no  errors  are  detected,  file  GCP  STORAGE  is  updated  using  a  direct  access  write  to 
GCP  STORAGE. 

ELMO  operates  in  a  manner  similar  to  SEMO  exeept  that  only  one  element  of  the  section 
is  modified  It  has  the  same  safeguards  as  SEMO  to  insure  that  the  proper  branch,  section, 
and  element  are  to  be  modified. 

ELMO  contains  a  list  of  element  types  in  logical  if  statement  and  these  are  scanned  and 
the  appropriate  location  of  the  STOR  array  is  updated  when  a  match  is  found.  ELMO  also  up¬ 
dates  GCP  STORAGE  as  SEMO  does. 
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I).  Diagnostic  Commands 

The  1, 1ST  command  calls  subroutine  QL1ST  which  prints  the  circuit  on  the  typewriter  console. 
Do  loops  cycle  through  ST  OH  (N,  M,  1,  2,  1)  (identification  code)  of  each  section  and  determine 
the  section  type  by  comparing  this  STOR  array  location  to  a  list  of  subroutine  identification  codes 
using  logical  ifs  When  a  match  is  found  the  if  action  is  a  write  statement  which  prints  the  section 
type. 

At  the  start  of  each  loop  the  do  loop  indices  are  printed  as  the  branch  and  section  numbers 
(N,  M).  The  indices  are  integer  variables  and  may  be  fixed  by  the  user  to  display  the  contents 
of  only  one  section.  This  is  done  in  two  ways. 

If  the  user  issues  the  commands  LIST  ALL,  ALL  is  read  into  the  integer  variable  ELEM1. 
This  variable  is  used  in  QL1ST  using  a  logical  if  to  set  the  do  loop  indices  so  that  all  sections 
and  branches  are  cycled  through. 

If  the  user  does  not  specify  ALL,  then  a  write  statement  asks  for  a  branch  or  a  branch  and 
section.  The  branch  and  section  is  read  into  two  integer  variables  and  these  are  used  as  the  do 
loop  indices.  If  just  a  branch  is  given,  the  do  loop  indices  are  fixed  so  that  the  whole  branch 
is  cycled 

A  nested  do  loop  is  used  to  display  the  values  of  the  section  elements  (i.e.,  the  third  dimen¬ 
sion  of  the  STOR  array).  Here  again,  if  zeroes  are  found,  the  do  loop  is  terminated  to  prevent 
printing  of  zeroes  for  sections  which  have  only  one  element. 

The  MATR  command  is  similar  to  the  LIST  command  except  that  the  ABCD  array  is  used 
instead  of  the  STOR  array,  and  the  printing  is  different.  It  is  handled  in  the  same  way  as  LIST 
as  far  as  the  do  loops  are  concerned. 

OFFL  sets  an  integer  variable  NR1TF  to  8.  All  write  statements  W'hich  w'rite  measured 
data  use  NR1TE  as  a  data  reference  set  number.  When  NRITE  equals  8,  all  writes  are  direct 
access,  and  write  on  a  file  named  GCP  DATA.  Thus  after  terminating  the  program  an  OFFLINE 
PRINT  GCP  DATA  must  be  issued. 

ONLI  sets  NRITE  to  6,  and  subsequent  writes  are  on  the  typewriter. 

The  final  command  to  GCP  is  QUIT.  This  command  terminates  the  execution  of  the  program. 
This  is  accomplished  by  scanning  the  list  of  commands  using  logical  ifs.  When  a  match  is  found, 
the  if  action  is  STOP. 

E.  Measurements 

The  MEAS  command  initiates  the  following  complex  series  of  actions. 

The  logical  ifs  are  scanned  for  a  match  to  MEAS.  When  it  is  found,  control  is  passed  to  a 
section  of  GCP  which  first  initializes  the  three  measurement  variables  to  blanks  (up  to  3  measure¬ 
ments  can  be  called). 

GCP  reads  the  measurement  into  3  integer  variables,  MEAS1,  MEAS2,  and  MEAS3,  and 
the  branch  and  section  into  Nl,  Ml  respectively.  The  variables  are  cleared  to  blanks  to  prevent 
any  previously  called  measurement  from  being  called  again. 

Section  and  branch  checks  are  performed  testing  STOR  (N,  M,  1,  2,  1)  for  zero  to  insure  a 
valid  measuring  point  (i.e.,  no  measurement  can  be  made  of  a  section  which  has  not  been  formed). 
An  error  message  is  printed  if  an  error  is  detected  and  control  returned  to  the  user  for  correc¬ 
tive  action.  If  no  error  is  detected,  the  subroutine  MEASUR,  which  now  acts  as  the  measure¬ 
ment  exec,  is  called,  passing  arguments  MEAS1,  MEAS2,  and  MEAS3. 
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Measurement  exec  first  checks  MEAS1,  MEAS2,  and  MEAS3  to  be  sure  they  are  not  all  blank 
using  logical  ifs.  If  they  arc  all  blank,  an  error  message  is  printed  and  control  is  returned  to 
command  exec  and  then  to  the  user  for  corrective  action. 

The  frequency  range  is  asked  for  with  a  write  statement.  \  NAME  LIST  HEAD  reads  the 
frequency  data  and  is  checked  to  insure  proper  form  (i.e.,  start  frequency  smaller  than  stop 
frequency  and  increment  not  equal  to  zero).  If  errors  are  detected,  appropriate  error  messages 
are  printed  and  control  is  passed  back  to  the  user  for  corrective  action. 

Integer  variables  NCHECK  and  NUMBER  are  set  to  zero.  NCHECK  is  used  to  suppress 
titles  (zero  prints  titles;  1  or  2  does  not  print  titles).  The  NUMBER  variable  is  a  counter 

Logical  ifs  are  used  to  determine  how  many  measurements  were  called.  This  is  accom¬ 
plished  by  cheeking  to  see  if  MEAS2  or  MEAS3  arc  blank.  If  more  than  two  measurements  were 
called  (MEAS2  not  blank)  then  the  subroutine  TITLE  is  called,  passing  the  arguments  MEAS1 , 
MEAS2,  and  ME  A  S3. 

This  subroutine  checks  the  list  of  measurements  against  the  passed  arguments  using  logical 
ifs,  and  writes  the  appropriate  titles  for  the  columns  on  the  typewriter 

Control  is  passed  back  to  measurement  exec  and  the  subroutine  QCALCZ  is  called.  This 
subroutine  controls  the  generation  of  the  transfer  matrices.  Using  two  nested  do  loops  it  cycles 
through  the  first  and  second  dimension  of  the  STOR  array  (these  locations  contain  branch  and 
section  numbers).  Each  time  it  goes  through  the  do  loop  it  calls  ZSUB  (an  entry  of  subroutine 
SUB),  which  checks  the  section  routine  identification  codes  using  logical  ifs.  When  a  match  is 
found,  the  "  Z"  entry  of  the  section  subroutine  is  called  and  this  in  turn  generates  the  transfer 
matrix  of  that  section  and  stores  it  in  the  ABCD  array. 

Control  returns  to  ZSUB  w-hich  returns  to  ZCALCZ,  and  the  do  loops  cycle  again  for  the 
next  section.  When  all  the  sections  to  be  measured  have  been  cycled  through  and  have  transfer 
matrices,  the  do  loop  is  exited.  This  is  accomplished  by  checking  STOR  (N,  M,  1,  2,  1)  for  zero 
as  before 

All  transfer  matrices  have  now  been  generated,  and  next  must  be  multiplied.  The  subroutine 
MATMUL  is  called  to  multiply  each  transfer  matrix  by  the  product  of  the  previous  matrices. 

MATMUL  is  simply  a  small  subroutine  for  multiplying  complex  matrices  by  the  product  of 
the  previous  matrices.  This  is  accomplished  with  a  do  loop  and  temporary  storage.  The  limits 
of  the  do  loop  are  set  by  QCALCZ.  The  upper  limit  Ml  being  carried  through  from  measurement 
exec  is  in  labeled  common  FREQ.  When  the  do  loop  is  satisfied,  control  returns  to  QCALCZ 
and  the  final  matrix  elements  are  stored  in  the  complex  variables  A,  B,  C,  and  D.  The  limits 
of  the  do  loop  (in  QCALCZ)  are  reset  to  provide  transfer  matrices  and  calculations  for  the  sec¬ 
tions  to  the  left  of  the  measuring  point  back  to  the  generator.  This  transfer  matrix  is  required 
for  some  measurements.  The  elements  of  this  matrix  are  stored  in  the  complex  variables 
AG,  BG,  CG,  and  DG. 

The  scattering  matrix  elements  are  also  generated  and  stored  in  SI  1 ,  SI  2,  S21,  and  S22. 
Control  then  passes  back  to  the  measurement  exec. 

Depending  on  the  number  of  measurements  called,  appropriate  action  is  performed  to  cor¬ 
relate  the  titles  and  columns  to  the  data  to  be  printed.  The  subroutine  SUBMEA  is  now  called. 
This  subroutine  scans  a  list  of  measurements  using  logical  if  to  find  the  measurement  that  has 
to  be  called.  When  a  match  is  found,  a  call  to  the  appropriate  entry  point  in  QMEAS  is  made. 
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QMEAS  is  a  subroutine  consisting  mainly  of  entry  points,  each  one  of  which  performs  cal¬ 
culations  to  obtain  the  measurement  requested.  The  result  of  the  calculation  is  stored  in  an 
array  ARR,  where  even  numbered  subscripts  contain  the  real  part  and  odd  numbered  subscripts 
contain  the  imaginary  part,  or  blank  if  the  result  is  all  real.  This  arrangement  is  necessary 
for  printing  purposes. 

If  only  one  measurement  is  requested,  QMEAS  will  print  the  results  before  passing  control 
to  the  measurement  exec.  If  more  than  one  measurement  is  called,  control  is  passed  back  to 
MEASUR  without  printing.  The  variable  NCHECK  is  used  to  test  whether  QMEAS  does  any 
writing. 

If  more  than  one  measurement  is  requested,  the  measurement  exee  will  print  out  the  results 
of  the  calculation.  The  results  of  the  calculations  are  in  the  labeled  common  MEA. 

After  printing  the  results  the  frequency  is  incremented  and  the  whole  process  is  repeated 
for  the  next  frequency.  When  the  stop  frequency  is  reached,  control  is  passed  back  to  the 
command  exee  and  the  user. 
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APPENDIX 

SUBROUTINES 


QSER 


Program  identification  code  —  2 
Other  subroutines  called  —  None 


Type  statements  needed 


COMPLEX  ABCD  (50,50,4),  RESULT  (4) 
REAL  LEN 


Common  statements  needed  —  COMMON/lNDEX/M,  N,  ABCD,  STOR  (10,  50,  4,  2,  2), 

W,  RESULT,  ELEM1 
COMMON/ZMEAS/L,  LI 


Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 


ABCD  Matrix 

SERI,  RESI 
A  1 

B  R  R  is  in  ohms 
C  0 
D  =  1 
SERI,  CAPA 
A  -  1 

1 

B  -  — 77  C  is  in  pF 

JwC  F 

C  -  0 

D  =  1 

SERI,  INDU 

A  1 

B  =  jo;  L  L  is  in  nil 
C  0 
D  =  1 
SERI,  TEMS 
A  =  1 

B  -  jZQ  sin  *9/cos  *0- 


C  -  0 
D  -  1 

where  *0-  =  uol  \[7[L/ 1 1.80271  x  10^ 
1  =  length  in  inches 
ep  =  dielectrie  constant 
SERI,  TEML 

A  1 


B  =  -  jZ  eotan  -0* 
J  o 


C  -  0 
D  =  1 

where  *0-  is  as  above 
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SERI,  PRLC 
A  =  1 

B  -  l/R  +  j (co C  —  1/ojL)  R,  L,  C  as  above 
C  0 
D  1 
SERI,  SRLC 
A  1 


B  R  +  j(wL  +  1/uC)  R,  L,  C  as  above 
C  0 
D  =  1 
SERI,  PRES 
A  -  1 


}  _  1 

B  =  BP  [1/Q  +  j(F/F0-F0'F)] 

C  0 
D  1 


F  =  frequency 

F0  =  resonant  frequency 

BP  =  susceptance 


SERI,  SRES 


A  1 

B  =  X  [1/Q  +  j(F/F  -  F  /F)]  F,  Fq,  BP  as  above 


C  0 
D  1 


QSER  is  a  subroutine  which  handles  all  series  elements.  The  element  type  is  passed  to 
QSER  stored  in  the  integer  variable  ELEM1  through  the  labeled  common  INDEX.  Using  a  list 
of  logical  ifs  an  appropriate  branch  is  made  to  a  part  of  the  subroutine  which  asks  for  the  in¬ 
formation  pertaining  to  the  element  type  and  stores  it  in  the  STOR  array. 

The  storage  location  STOR  (N,  M,  1,  1,  2)  is  eoded  for  use  when  generating  the  transfer 
matrix. 


II.  QSHU 


Program  identification  code 
Other  subroutines  called 
Type  statements  needed 

Common  statements  needed 


3 

None 

COMPLEX  A  BCD  (10,50,4),  RESULT  (4). 

REAL  LEN 

COMMON/INDEX/M,  N,  ABCI),  STOR  (10,50,4,2,2), 
W,  RESULT,  ELL  Ml 
COMMON/ZMEAS/L,  LI 


Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCD  Matrix 


SHUN,  RESI 
A  =  1 


B  =  0 

C  R  R  is  in  ohms 
D  1 
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SHUN,  CAPA 
A  =  1 
B  0 

C  ju;C  C  is  in  pF 
D  1 
SHUN,  INDU 
A  =  1 
B  0 

1 

C  - — t~  L  is  in  nH 

jcjL 

D  1 

SHUN,  TEMS 

A  -  1 

B  -  0 

COS  *9- 

c  ~  iZ  sin  -0 
J  o 

D  =  1 

**  Q 

11.80271  X  10y 

where  1  =  length  in  inches 

ep  dielectric  constant 

SHUN,  TEML 

A  =  1 

B  -  0 

r  _  tan  *©■ 

lZ 
J  o 

D  =  1 

where  #  and  e m-  are  as  above 
SHUN,  PRLC 
A  1 
B  -  0 

C  -  R  +  i(wC  -  -V)  R,  C,  L  as  above 

CaJ  Li 

D  1 
SHUN,  SRLC 
A  -  1 
B  =  0 

C  4  +  j(wL  +  — U)  R,  C,  L  as  above 

K  CaJ  L, 

D  =  1 
SHUN,  PRES 
A  =  1 
B  =  0 


D  =  1 


29 


SHUN,  SUES 
A  1 
B  0 


r _ 7 _ 1 _ 

X[l/Q+  j(F/F0-F0/F)] 

D  1 

QSHU  is  a  subroutine  whieh  handles  all  shunt  elements,  and  is  identical  to  QSER  exeept  for  the 
generation  of  the  transfer  matriees. 

III.  QTEML 

Program  identifieation  code  —  4 
Other  subroutines  called  —  None 

Type  statements  needed  —  COMPLEX  ABCD  (10,50,4) 

HEAL  LEN 

Common  statements  needed  -  COMMON/lNDEX/M,  N,  ABCD,  STOR  (10,  50,  4,  2,  2),  W 

COMMON/ZMEAS/L,  LI 

Returns  -  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCI)  Matrix 

A  =  eos  -0- 

B  =  iZ  sin  -0- 
J  o 

C  - - - 

iZ  sin  -0- 
J  o 

I)  =  eos  O 

where  -0-  =  o;l  \fe\±/i  1 .8  027 1  x  10^ 

1  =  length  in  inehes 
€\x  dielectric  eonstant 

This  subroutine  simulates  a  length  of  transmission  line.  Zq,  length  and  dielectrie  eonstant 
are  read  into  variables  Z0,  LEN,  and  EPMU,  then  stored  in  the  STOR  array. 

IV.  QSHORT 

Program  identification  code  —  5 
Other  subroutines  ealled  —  None 

Type  statements  needed  —  COMPLEX  ABCD  (10,50,4) 

Common  statements  needed  —  COMMON/lNDEX/M,  N,  ABCD,  STOR  (10,  50,4,2,2),  W 

COMMON/ZMEAS/L,  LI 

Returns  -  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCD  Matrix 

A  -  0 
B  1 
C  1 
I)  =  0 

This  subroutine  simulates  a  short  cireuit  at  a  braneh  end.  The  user  docs  not  supply  any 
information  to  this  subroutine. 
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V.  ©JUNC 


Program  identification  code  —  6 


Other  subroutines  called 
Type  statements  needed 


-  ZZSUB,  Z 1 S  U  B ,  MATMUL 


-  COMPLEX  ABCD  (10,50,4),  RESULT  (4), 

ZIN,  YIN 
INTEGER  SER,  SHU,  CONN,  QUIT, 

BLK,  ELEM1,  TYPE 

Common  statements  needed  —  COMMON/lNDEX/M,  N,  ABCD,  STOR  (10,  50, -4,  2,  2), 

W,  RESULT,  ELEM1 
COMMON/ZMEAS/L,  LI 

CO  M  MON/ FREQ/ST  A  F,  STOF,  DELF,  Ml,  Nl,  N2, 
N3,  M2,  M3 

Returns  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCD  Matrix 

SHUNT  JUNCTION 
A  =  1 
B  0 


C  =  Y. 


in 

D  1 

SERIES  JUNCTION 
A  1 


Y^  is  C/A  of  the  produet  of  transfer  matrices  of  the 
specified  branch. 


B  -  Z. 


in 


Zin  is  A/B  of  the  produet  of  transfer  matrices  of  the 
specified  braneh. 


C  -  0 
D  =  1 

This  subroutine  is  used  to  attach  SHUNT  or  SERIES  branches  to  the  main  branch,  or  to 
another  branch. 

The  subroutine  has  three  sections.  The  first  section  deals  with  the  mechanics  of  requesting 
new  branches  and  entering  a  branch  at  the  next  available  section. 

The  seeond  and  third  sections  deal  with  the  attaching  of  branches  to  other  branches.  The 
third  section  is  actually  a  separate  subroutine  whieh  is  required  because  multiple  calls  to  the 
same  subroutine  are  not  possible. 

Branch  and  section  numbers  are  stored  in  STOR  (N,  M,  1,  1,  1)  and  STOR  (N,  M,  2,  1,  1)  re¬ 
spectively.  Logical  ifs  test  the  variable  type  (containing  the  alphanumeric  SHUN  or  SERI)  to 
determine  the  connection  of  the  branch  (SHUNT  or  SERIES).  The  constant  1  is  stored  in  STOR 
(N,  M,  3,  1,  1)  if  a  SHUNT  connection  is  desired,  and  2  is  stored  if  a  SERIES  connection  is 
desired. 

To  generate  the  transfer  matrix  for  this  section,  it  is  necessary  to  call  the  subroutine 
ZZSUB  which  performs  the  same  function  as  ZSUB  described  in  Sec.  IV-E.  This  subroutine 
generates  the  transfer  matrices  for  the  seetion  of  the  braneh  to  be  attached.  When  all  the  trans¬ 
fer  matrices  have  been  formed,  MATMUL  is  ealled  to  obtain  the  product  of  the  matrices.  Then 
the  appropriate  transfer  matrix  is  generated  for  the  junction. 

To  attach  a  branch  to  a  branch,  subroutines  ZZJUNC  and  Z1SUB  are  used  as  above  for  a 
single  branch  connection. 


31 


VI.  QEMPTY 

Program  identification  code  —  7 
Other  subroutines  called  —  None 

Type  statements  needed  —  COMPLEX  ABCD  (10,  50,4) 

Common  statements  needed  -  COMMON/lNDEX/M,  N,  ABCD,  STOR  (10,  50,4,  2,  2) 

COMMON/ZMEAS/L,  LI 

Returns  -  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCD  Matrix 

A  1 
B  -  0 
C  -  0 
D  -  1 

This  subroutine  consists  of  just  an  identity  matrix  which  simulates  an  empty  section. 

VII.  QISOL 

Program  identification  code  —  8 

Other  subroutines  called  —  None 

Type  statements  needed  —  REAL  LF,  LR 

COMPLEX  ABCD  (10,50,4) 

Common  statements  needed  -  COMMON/lNDEX/M,  N,  ABCD,  STOR  (10,  50,4,  2,  2) 

COMMON/ZMEAS/L,  LI 

Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCD  Matrix 

A  =  (x  +  y) 

B  =  Zq  (x-y) 

C  -  - - - 

zo 

D  =  (x  +  y) 

10LF/2° 

where  x  —  LF  is  forward  loss  (dB) 


y  ~  - t  p  /7a  LR  is  reverse  loss  (dB) 

2  (10I"K  Z0) 

This  subroutine  simulates  an  isolator.  The  forward  and  reverse  loss  in  dB  is  read  into 
LF  and  LR,  then  stored  into  STOR  (N,  M,  2,  1,  1)  and  STOR  (N,  M,  3,  1,  1)  as  10.**  (LF/20)  and 
10.**  (LR/20),  respectively.  ZO  is  stored  in  STOR  (N,  M,  1,  1,  1). 


QATTN 

Program  identification  code 

-  9 

Other  subroutines  called 

—  None 

Type  statement  needed 

-  REAL  L 

COMPLEX  ABCD  (10,  50,  4) 

Common  statements  needed 

-  COMMON/INDEX/M,  N,  ABCD,  STOR  (10,  50,  4,  2,  2) 
COMMON/ZMEAS/L,  LI 
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Returns 


-  Statement  No.  1  of  SUB 
Statement  No.  5  of  QSEMO 

A  BCD  Matrix 

A  (x  +  v) 

B  Zq  (x  -  v) 

C  - - - 

zo (x  - 

D  (x  +  y) 

where  x  -  i  (10^  L  -  Attenuation  (dB) 

2  [10(L/20)] 

This  subroutine  simulates  a  matched,  fixed,  attenuator.  The  attenuation  is  read  into  L 
then  stored  into  STOR  (N,  M,  2,  1,  1)  as  10.**  (L/20).  is  stored  in  STOR  (N,  M,  1,  1,  1). 

IX.  QXFMR 

Program  identification  code  —  10 
Other  subroutines  ealled  —  None 

Type  statements  needed  —  COMPLEX  ABCD  (10,  50,4) 

Common  statements  needed  -  COMMON/lNDEX/M,  N,  ABCD,  STOR  (10,  50,  4,  2,  2) 

COMMON/ZMEAS/L,  LI 

Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCD  Matrix 

A  =  N 
B  0 
C  0 

D  =  il 

This  subroutine  simulates  an  ideal  transformer.  The  turns  ratio  is  read  into  TURNS  then 
stored  in  the  STOR  array. 


X.  QINVE 


Program  identification  code 
Other  subroutines  ealled 
Type  statements  needed 
Common  statements  needed 


11 

None 

COMPLEX  ABCD  (10,50,4) 

COMMON/INDEX/M,  N,  ABCD,  STOR  (10,50,4,2,2) 
COMMON/ZMEAS/L,  LI 


Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 


ABCD  Matrix 


A  0 
B  Z 

o 


C 


1 


jZ 


o 


D  0 
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This  subroutine  simulates  an  impedance  inverter.  ZQ  is  read  into  ZO  then  stored  in  the 
STOR  array. 


XL  QCROS 

Program  identification  code  —  12 
Other  subroutines  called  -  None 

Type  statements  needed  —  COMPLEX  ABCD  (10,50,4) 

Common  statements  needed  -  COMMON/lNDEX/M,  N,  ABCD,  STOR  (10,  50,  4,  2,  2) 

COMMON/ZMEAS/L,  LI 

Returns  -  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCD  Matrix 

A  =  —  1 
B  -  0 
C  -  0 
D  =  -1 

This  subroutine  simulates  a  change  in  polarity.  No  information  is  required  from  the  user. 


XII.  QRADL 


Program  identification  code 
Other  subroutines  called 
Type  statements  needed 

Common  statements  needed 


13 

BESJ,  BESY 

REAL  J1X,  JOX,  J1Y,  JOY 
COMPLEX  ABCD  (10,50,4) 

COMMON/INDEX/M,  N,  ABCD,  STOR  (10,50,4,2,2),  W 
COMMON/ZMEAS/L,  LI 


Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 


ABCD  Matrix 


A  -  [J  1(Y)  Y 0(X)  -  Y1(X)  J0(X)J 


B  =  LX  [ J0( Y)  Y0(X)  -  Y0(Y)  J0(X)] 


C  -  [J1(Y)  Y1(X)  —  Y1(Y)  J1(X)J 

D  -  [Y0( Y)  J1(X)  —  J0(Y)  Y1(X)J 

ZRo 


where  X  = 


Z 


\[T\ I  c o  R 
11.8  x  109 
60B 


°  R  \TejI 
o 

B  is  the  spacing 

R  is  radius  at  the  generator  end 

R  is  radius  at  the  load  end 
o 

C(jl  is  the  dielectric  constant 

Jl,  JO,  Yl,  Y0  are  BESSEL  functions 
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This  subroutine  simulates  a  radial  line.  The  diameter  of  the  load  and  generator  ends  are 
read  into  Rq  and  R,  respectively.  The  spacing  and  dielectric  constant  are  read  into  B  and 
EPMU,  respectively.  STOR  (N,  M,  1,1,1)  contains  R/2.,  STOR  (N,  M,  2,  1,  1)  contains  EPMU**. 5, 
SI  OR  (N,  M,  3 ,  1 ,  1 )  contains  B,  and  STOR  (N,  M,  4,  1,  1)  contains  R^  Z. 

The  entry  ZRADL  calculates  X  and  V  as  above  and  calls  the  various  BESSEL  subroutines. 

The  BESSEL  subroutines  return  an  error  flag  and  this  is  printed  if  an  error  occurs.  The 
following  table  indicates  the  errors: 

Error  Code  Error 

Z  X  is  negative 

3  Accuracy  not  obtained  (set  at  one  percent); 

or  returned  argument  >10 ' 


XIII.  QDUMM 

Program  identification  code 


Other  subroutines  called 
Type  statements  needed 
Common  statements  needed 


QDUMMA  -  14 
QDUMB  -  15 
QBUMC  -  16 
ODUMD  -  17 

None 

COMPLEX  (10,  50,  4) 

COMMON/INDEX/M,  N,  ABCD,  STOR  (10,50,4,3,3),  W 
COMMON/ZMEAS/L,  LI 


Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCD  Matrix 

As  required  by  user. 


This  subroutine  consists  of  four  entries  which  the  user  may  use  to  write  his  own  programs. 


XIV.  Ql  MAG 

Program  identification  code  —  18 
Other  subroutines  called  —  None 

Type  statements  needed  —  COMPLEX  ABCD  (10,  50,4),  TEMP  (4),  A,  B,  C,  D 

Common  statements  needed  -  COMMON/lNDEX/M,  N,  ABCD,  STOR  (10,  50,4,2,2) 

COMMON/ZMEAS/L,  LI 

Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCD  Matrix 

Product  of  several  matrices. 

This  subroutine  simulates  a  mirror  image  of  specified  sections.  Start  and  end  branch  and 
section  numbers  are  read  into  NS,  MS,  NE,  and  ME,  then  stored  in  the  STOR  array. 

The  transfer  matrix  is  generated  as  in  MATMUL,  writh  the  product  stored  in  the  ABCD  array. 


XV.  QTEWG 

QTEWG  is  an  entry  point  in  QTEML. 

Program  identification  code  —  4 
Other  subroutines  called  —  None 
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Returns 


-  Statement  No.  1  of  SUB 
Statement  No.  5  of  QSEMO 

ABCD  Matrix 

A  cos  0 

B  =  iZ  sin  0 
J  o 

C  - - - 

jZ  sin  0 
o 

D  cos  0 
where 

e  2tt  7 (FREQ1Z  -  FC2)  (LEN)  n/T^T 
11.80271  x  lo9 

FREQ1  § 

FC  11.80271  X  109/[2(\UDE)  \/l>] 
WIDE  =  wide  dimension  of  waveguide 

z  _ _ 120a;  (NARO) _ 

WIDE  477  J( FREQ12  -  FC2) 
NARO  =  Narrow  dimension  of  waveguide 


XVI.  QWLOD 

QWLOD  is  an  entry  point  in  QTEML. 

Program  identification  code  —  20 

Other  subroutines  called  —  None 

Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 

ABCD  Matrix 

A  =  1 
B  0 
C  -  Y 

o 

D  1 
where: 

(WIDE)  47^  /  (^-)2  -  FC2 
Y  -  _ V  _ 

o  120oj  (NARO) 

NARO  =  narrow  dimension  of  waveguide 
WIDE  =  wide  dimension  of  waveguide 


XVII.  QWISO 

QWISO  is  an  entry  point  in  QISOL. 

Program  identification  code  —  21 

Other  subroutines  called  -  None 

Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 
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ABCD  Matrix 


A  i/2(FL  +  j^-) 


2(FL"  RL> 


Y 


C  = 


2<fl-rl> 

D  =  1/2(FL  +  ££) 


where 


FL  -  10 

( 

RL  =  10 
Z 


, forward  loss, 
1  20  } 


, reverse  loss, 

1  20  } 

1  20co  (NARO) 


(wide)  kRJL  J(%)Z  -  FC" 


o  Z 


FC  = 


11.80271  X  10' 


2<WIDE)  \T iT 

WIDE  =  wide  dimension  of  waveguide 
NARO  narrow  dimension  of  waveguide 


XVIII.  QPARB 


Program  identification  eode  —  22 
Other  subroutines  called  —  None 


Type  statements  needed 
Common  statements  needed 


COMPLEX  ABCD  (10,  50,4) 

INTEGER  ELEM1 

COMMON/INDEX/M,  N,  ABCD,  STOR  (10,  50,  4,  2,  2), 
W,  RESULT,  ELEM1 
COMMON/ZMEAS/L,  LI 


Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 


ABCD  Matrix 


A  1 


B  =  0 
C  0 
D  1 


This  subroutine  signals  the  start  of  a  parallel  branch. 

The  command  is  PARB  p  where  p  is  a  branch  number.  The  branch  number  is  read  into 
the  integer  variable  ELEM1  in  an  A4  format.  It  is  then  written  on  data  reference  set  number  1 
in  an  A 4  format  and  read  baek  as  an  12  format.  This  is  required  to  convert  it  from  an  A  format 
to  an  I  format. 

The  braneh  number  is  stored  in  STOR  (N,  M,  1,  1,  1)  to  be  used  later  in  the  QCOMB  subroutine. 
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With  logical  ifs  the  branch  number  is  checked  to  make  sure  that  a  branch  number  was  given, 
and  that  it  is  between  1  and  10. 

The  A  BCD  matrix  is  set  to  the  unit  matrix  because  the  program  computes  the  products  of 
the  ABCD  matrix  of  the  section  when  a  measurement  is  called  for. 


XIX.  QCOMB 


Program  identification  code 
Other  subroutines  called 

Type  statements  needed 


Common  statements  needed 


MATMUL 

Z1SUB 

COMPLEX  ABCD  (10,  50,4),  RESULT  (4), 

TEMP  (4),  DELI,  DEL,  ZDEL,  Zll,  ZZ11, 

Z12,  ZZ12,  Z21 ,  ZZ21,  Z22,  ZZ22 

COMMON/INDEX/M,  N,  ABCD,  STOR  (10,50,4,2,2) 
W,  RESULT 

COMMON/ZMEAS/L,  LI 

COMMON/FREQ/STAF,  STOF,  DELF,  Ml,  Nl,  N2, 
N3,  M2,  M3 


Returns  —  Statement  No.  1  of  SUB 

Statement  No.  5  of  QSEMO 


ABCD  Matrix  (both  paths  not  pure  shunt,  see  below) 

*  B'A"  +  B"  A' 

B’  4  B" 


B 


B'B" 
Bf  +  B" 


C  =  C'  +  CM  4  AM 
D  =  1  - 


Dn 


rv  D’  —  Dn 

)  +  A»(p,  x  g„) 


B1  4  B"'  '  “  'B’  4  B 

Bf  ( 1  -  D")  4  BM  (1  -  Df) 

Br  4  B" 


where: 

A1,  B’,  C1,  D1  are  the  elements  of  the  matrix  formed  by  multiplying  the 
transfer  matrices  of  one  parallel  branch. 

AM,  B",  C",  D"  are  the  elements  of  the  matrix  formed  by  multiplying  the 
transfer  matrices  of  the  other  parallel  branch. 

ABCD  Matrix  (both  branches  are  pure  shunt) 

A  =  1 


B  -  0 


c  -  C’  4  C" 

D  -  1 

where  C’,  C"  are  as  above. 


QCOMB  performs  the  function  of  creating  transfer  matrices  for  both  branches,  finding  their 
products,  and  combining  them  into  a  single  transfer  matrix. 

The  subroutine  first  saves  the  integer  variables  Ml,  M2,  M3,  Nl,  N2,  and  N3,  because  the 
original  values  are  destroyed. 

Using  a  do  loop  the  main  branch  is  searched  to  determine  the  number  of  parallel  sections 
it  contains.  This  is  done  by  checking  STOR  (L,  LI,  1,  2,  1)  for  22,  which  is  the  program  identifi¬ 
cation  code  for  PARB  indicating  the  start  of  the  parallel  branch. 

After  this  is  determined  the  subroutine  MATMUL  is  called  to  multiply  the  transfer  matrices 
of  these  sections,  and  get  product  matrix.  The  integer  variables  M2,  M3  are  set  to  the  proper 
values  so  that  MATMUL  multiples  just  the  sections  that  are  in  parallel  with  the  other  branch. 
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The  other  parallel  branch  is  then  searched  to  determine  the  number  of  sections  by  checking 
STOR  (10,  50,  4,  2,  Z)  for  zeroes  in  a  do  loop. 

The  subroutine  Z1SUB  is  called  to  generate  the  transfer  matrices  for  these  sections,  and 
MATMUL  is  called  to  find  their  product. 

At  this  point  the  product  matrix  of  both  branches  is  checked  to  determine  if  both  branches 
are  pure  shunt.  If  both  branches  are  pure  shunt  the  B  element  is  zero;  checking  the  ABCD 
matrix  above,  it  can  be  seen  that  a  zero  divide  exception  will  occur  if  this  matrix  is  generated. 
Using  a  logical  if,  the  B  element  of  each  branch  is  checked  for  zero,  and  if  they  are  not  both 
zero,  the  variables  Ml,  Nl,  M2,  M3,  N2,  N3  are  returned  to  their  original  values  and  the 
matrices  arc  combined  as  in  the  first  A  BCD  matrix. 

If  the  B  elements  are  both  zero  the  matrices  are  combined  as  in  the  second  ABCD  matrix 
above. 

Since  the  matrices  of  the  main  branch  are  multiplied  to  form  a  product  matrix  from  which 
the  measurements  are  made,  it  is  necessary  to  change  the  sections  of  the  main  branch  which 
are  in  parallel  with  the  second  branch  to  unit  matrices.  This  is  done  using  a  do  loop. 
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